Exact decoherence to pointer states in free open quantum systems is universal 
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In this letter it is shown that exact decoherence to minimal uncertainty Gaussian pointer states is generic 
for free quantum particles coupled to a heat bath. More specifically, the paper is concerned with damped free 
particles linearly coupled under product initial conditions to a heat bath at arbitrary temperature, with arbitrary 
coupling strength and spectral densities covering the Ohmic, subohmic, and supraohmic regime. Then it is true 
that there exists a time t c such that for times t > t c the state can always be exactly represented as a mixture 
(convex combination) of particular minimal uncertainty Gaussian states, regardless of and independent from the 
initial state. This exact 'localisation' is hence not a feature specific to high temperatures and weak damping 
limit, but is rather a generic property of damped free particles. 
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There is long tradition of approaching the questions of how 
and to what extent classical properties of quantum systems 
emerge dynamically due to the unavoidable coupling to their 
environment. Essentially any quantum system interacts to 
some extent with other external degrees of freedom, which in 
turn may be said to monitor certain properties of the quantum 
systems 

QllSllEEHSIlIli. This y ields decoher- 
ence, which results in a loss of purity of initially pure states of 
a distinguished quantum system coupled to an environment. 
Not all initial quantum states in such a dilation are yet equally 
'fragile' to this interaction: there is a small set of initial states 
that is often relatively robust with respect to this interaction. 
The term pointer states has been coined for such states, ow- 
ing the name to models for quantum measurement where the 
pointer basis is essentially determined by the interaction of the 
apparatus with the external degrees of freedom 

For harmonic and free quantum systems linearly coupled 
to a heat bath consisting of harmonic systems this general 
mechanism is very well-understood indeed. For example, if 
one prepares a single mode in a pure state in order to let it 
very weakly interact with an environment in the Gibbs state 
corresponding to a very high temperature, which one is the 
state that produces the least entropy over one cycle of the 
oscillator? In retrospect it hardly comes as a surprise that 
this is a coherent state, a minimal uncertainty Gaussian state 
0] . Most attention has probably been devoted to thoroughly 
understanding the dynamics of harmonic and free quantum 
systems in this limit of weak coupling and high temperatures 
GIBBS]- In this limit in particular, decoherence time scales 
have been identified 0] . But also exact quantum master equa- 
tions, generators of dynamical maps, have been derived and 
scrutinized in great detail [HI E B Hi- After all, the 
dynamics of open harmonic or free quantum systems can not 
be described other than being well-understood. What else is 
there to ask for? 

A question that seems to have been overlooked so far, yet, 
is the following: To what extent is exact decoherence in free 
quantum systems to pointer states actually generic? This 
question is most appealing in case of the free damped quan- 
tum particle fl2l uM . where there is no equilibrium Gibbs 
state. More specifically: is it true that starting from an ar- 



bitrary quantum state, after a fixed finite time t c (independent 
of the initial state), the state of the system is exactly indis- 
tinguishable from a mixture, a convex combination, of mini- 
mal uncertainty Gaussian states for all times t > t c 7 In this 
sense the free quantum system may be said to be in a situation 
that can operationally not be distinguished locally from the 
following situation: the particle is somewhere, in a minimal 
uncertainty Gaussian state, one simply does not know where 
in phase space. That this is the case seems fairly plausible 
for the case of high temperatures and weak damping. A sig- 
nificant first step in this direction has indeed been achieved 
very recently by Diosi and Kiefer in Ref. 01, showing that 
this intuition is indeed correct for the approximate generator 
for the dynamical map in the limit of negligible friction and at 
high temperatures. But is this a generic feature of free quan- 
tum systems that are linearly coupled to an environment in 
a dilation, and true not only for specific regimes, but for any 
coupling strength, any non-zero temperature, and Ohmic, sub- 
ohmic, as well as supraohmic damping? This is the question 
that will be addressed (and answered) in this paper. 

A free quantum system linearly coupled to a heat bath of os- 
cillators will be investigated, where the distinguished system 
is initially in an arbitrary (and potentially very 'non-classical' 
state), whereas the environment is prepared in the Gibbs state, 
which corresponds to an initial product state, such that the 
time evolution of the state of the free quantum system amounts 
to a completely positive dilation Hal . No assumptions will be 
made concerning the temperature of the environment and the 
strength of the coupling; for the class of non-vanishing spec- 
tral densities any C°° -function / : R + — ► R+ could be 
allowed for with 

lim I(uj)/uj p = C > (1) 

for some p € (0, 2). This will be referred to as Ohmic damp- 
ing when p = 1, otherwise as subohmic (for p < 1) or 
supraohmic (for p > 1). This is an already solved prob- 
lem in the sense that quantum master equations are known, 
and hence, the argument draws heavily from known results 
on generators of dynamical maps II ll 1 1711 . and from earlier 
results on the long-time behavior in quantum Brownian mo- 
tion (3- The starting point is the equation of motion of the 
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reduced density operators as derived in Ref. fiUl . in the inte- 
grated form as presented in the recent paper Ref. llTj7ll . Later, 
ideas will be used very similar to the ones in Ref. 0|. 

The equation of motion of the free particle is for the subse- 
quent purposes most conveniently be expressed in phase space 
in terms of the Wigner function W : R 2 x 1R + — > R Hal , 
which is for each t £ R + the Fourier transform of the char- 
acteristic function, dependent on £ = (^1,^2) €E R 2 , where 
£1 and £2 correspond to position and momentum coordinates 
in phase space, respectively. As a partial differential equation 
the Hu-Paz -Zhang equation fiUl reads lH9ll 

d t W(£, t) = -£2% t) + f2 2 (i)&% W(£, t) 
+2T(t)d €a faW{t,t)) 

+T(t)h(t)dl W{£, t) + r(t)/(t)% d i2 W{t t) , 

where the r, f,h,Q : R + — ► R are time-dependent coeffi- 
cients for which explicit expressions are known. The formal 
solution of this partial differential equation can be found for 
all system parameters The solution of the differen- 

tial equation with time-dependent coefficients as presented in 
Ref. ilVII is given by 



as G(0) = 0. Eq. (0, together with the subsequent specifica- 
tions forms the starting point of our analysis. 

Eq. can, using the transformation rule for multiple in- 
tegrals, be written in form of a product of a time dependent 
determinant and a convolution with a Gaussian as 



1 



1 



2 7 r|M(t)| 1 /2 
W(V(t)-^',0) 



\v(t)\ 

where the 2 x 2-matrix V is given by 



(2) 



V(t) 



G{t) G(t) 
G(t) G(t) 



1 



2w\M(t)\ 1 / 2 

where dots represent time derivatives. Here, G : R — > R is 
the Green's function, which is G(t) = for t < and is for 
t > the solution of the integral equation 



The Green's function G can not be evaluated in general in a 
closed form, the case of Ohmic damping being an exception, 
where the spectral density is for small frequencies linear in 
the frequencies. The Laplace transform G of G is related to 
the Laplace transform 7 of 7 as G(z) = [z 2 + zj(z))~ 1 . In 
order to specify the long time behavior of the Green's func- 
tions, it is sufficient to know the power law for the spectral 
density for small frequencies only. Using Eq. Q, one arrives 
forp 6 (0, 2) at lim^^ G(t)/f(t) = 1 (see also Ref. Q), 
where /(t) = sin(7rp/2)t p_1 /((T(p)). From the asymptotic 
behaviour of / as t — > 00 it can be seen after a few steps that 
lim^oo A(t)/A'(t) = 1, with 

A'(t)= f ds f ds'G(t-s)G(t-s')K(s-s'), 
J — 00 J 00 

This quantity in turn happens to be a quantity investigated in 



Ref. 



qua 

M 



G(t) + / ds7(i-s)G(s)=0,7(t) 
Jo 



aw cos(u)t), 



where it has been shown that 

\im(A'(z)(3z)/(2G(z)) = 1, 
2^0 



with initial conditions G(0) = 0, G(0) = 1, in terms of the 
so-called damping kernel. The 2 x 2-matrix 



M(t) 



A(t) C(t) 
C(t) B(t) 



which yields lim^oo A(t)/A"(t) = 1, with 

= 2sin(^/2) 
1 ; KT(p + 1) 



(3) 



has coefficients that have in Ref. fTvtl been expressed in terms 
of correlation functions. On using the function K : R + — ► 
R, 

1 f°° 

K{t) = — / duire[jiui + i0 + )] wcoth(/3w) cos(wt), 
t Jo 

with 7 = J °° dtjit)e lzt , the coefficients A(t), B(t), and C(t) 
can be expressed as 

A(t) = 

Bit) = 

C(t) = 



In order to find the long-time behaviour of the function C, we 
may use the fact that C{t) = 2A{t) for all t e [0, 00), which 
holds since G(0) = 0, and apply l'Hospital's rule to arrive at 
lim^oo C(t)/C"(t) = 1, with 



„ _ gpan(7rp/2) ! 



(4) 



ds 



ds 



s)Git - 


s')K(s - s 


s)G(t - 


s')Kis — s 


s)G(t - 


s')Kis - s 



To get the long term behaviour of B, we can again start with 
that lim^oc B{t)/B'{t) = 1, where 

B'(t) = I ds I ds'G(t - s)Git - s')K(s - s'). 

This, in turn, is nothing but the momentum uncertainty in the 
stationary setting, which is well-defined even in this free case 
(compare also Ref. PEHn . 

1 f°° 

- / duim{aiu) + i0 + )]uj 2 coth(w/3) =5 oo >0, (5) 
n Jo 
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with a(z) = —l/(z 2 + izj(z)). B^ is a (time-independent) 
positive real number. So we have determined the long-time 
behaviour of the entries of the symmetric 2 x 2-matrix M(t). 

Subsequently, a pointer state is taken to be a minimal uncer- 
tainty Gaussian state with particular second moments that re- 
flect a small uncertainty in the position canonical coordinate. 
The statements will be formulated in a language common in 
quantum optics and continuous-variable quantum information 
theory. The first moments are (d±, d%) = ((X), (P)), the sec- 
ond moments are collected in the covariance matrix 



2(0?) (OiOa + OaOi) 
<0i0 2 + O a Oi) 2(0 2 > 



where Oi = X - {X) and 2 = P - (P). 
moments for the pointer states are taken to be 



Boo 1 






Boa 



The second 



(6) 



This is a covariance matrix of a minimal uncertainty state, as 
Too | = 1. Note that in the weak damping limit, Boo becomes 
approximately 1231 -Boo = P~ 1 =T (note that K = 1 and k = 
1), so in the weak coupling and high temperature limit, the set 
of pointer states is a set of minimal uncertainty Gaussian states 
very narrow in position. The corresponding pure Gaussian 
state with first moments (di,d 2 ) = (£i , £2 ) will t> e denoted as 

This set of minimal uncertainty Gaussian states, which be- 
comes a set of states very narrow in position in the limit of 
weak coupling and high temperatures, will be regarded as the 
set of pointer states 12311 . It is an overcomplete set of states sat- 



isfying |("0dV>e)f 



-(£- 



=/2)(£-0 



The analogue 



of the standard s-ordered Wigner function of a state p may be 
defined as 

s E [— 1, 1], where a is the symplectic matrix embodying the 
canonical commutation relations, {.,.}+ denotes the anticom- 
mutator, and a = —\og(B 00 )/2 is the squeezing parameter 
corresponding to the pointer states (taken with respect to the 
standard unit quantum oscillator). The state can then be rep- 
resented as 12411 



P = / d 2 £WMOI^X^I, 



(7) 



whereas in turn W—i(£) = (i/j^\p\i/j^)/tt > for all £ £ 
R 2 . Then, the s-ordered functions are related to each other 
via convolutions (compare, e.g., Ref. O) 



-2(«-«')r- 1 («-«') T /(s'- s ) 



2tt 



(8) 



for s < s' , We are now in the situation that we can argue 
similarly to Ref. 0): the function Wq : R 2 x R+ — > R, 




-3-2.5-2-1.5-1-0.5 



lo 



FIG. 1: This figure shows logT c = lim-r^oo log t c , where for a 
given temperature the number t c > is the smallest number such 
that is satisfied for all t > t c for the case of strictly Ohmic damp- 
ing, as a function of log f , The stronger the damping, the faster is the 
'localization' process. 



W^,t) = Wo(V- 1 (t)£,0)/\V(t)\, is a legitimate Wigner 
function, as can be read off the definition of the Wigner func- 
tion. Then, Eq. (|2ji and Eq. © imply that 



Wi(f,t) 



x e 



d 2 r MM| M(t) _ roo/4 |-l/2 
27T 

-(?-«')(A/(t)-r 00 /4)- 1 (c-C) T /2 



But since 



rf 2 r MM 4e - 

2tt 



2(?-?')roo(«-«') T 



> 



for all £ e R 2 , Wi (£, t) > for all £ S R 2 if 

M{t) - Too/A > Too/A. 



(9) 



In turn, given the time dependence of the coefficients of M(t), 
demonstrated in Eqs. Q, @, and (0, there exists a finite t c > 

such that (|9} is valid for all t > t c . This time t c , in turn, is 
the time from which on W\ is strictly positive, and the state 
can certainly exactly be represented as a mixture of pointer 
states with second moments as in Eq. l|6). 

This is a generic result for arbitrary non-zero temperatures, 
arbitrary coupling strengths and all the spectral densities as in 
Eq. JO. For specific choices for the spectral density, bounds 
for the time t c can be found from which on the state can be 
represented as a mixture of pointer states. For Ohmic damp- 
ing in particular, the Green's function is given by G(t) = 
(1 - e" ct )/C, i.e., j(z) = ( > 00. The behaviour be- 
comes particularly transparent in the high temperature case. 
We then simply obtain lirriT^oo A{t)/T = 2{t - G(t))/(, 
lirriT^ C(t)/T = 2(1 - G(t))/(, and lim T ->oo B(t)/T = 

1 — e~ 2 ^*. Fig. ^ depicts T c = limy^oo t c , where t c is 
the smallest time for which © is satisfied for strictly Ohmic 
damping. 

To conclude, it has been shown that if one couples a free 
particle linearly to a heat bath prepared in the Gibbs state of 
some temperature, then, under very general conditions and 
without approximations, the state of the system becomes after 
some finite time exactly indistinguishable from an exact mix- 
ture of particular minimal uncertainty Gaussian pointer states. 
In this sense it can be said that exact decoherence to these lo- 
calized pointer states is generic, and not only a feature of a 
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limit that can be regarded as being classical. Locally, hence, 
we arrive at the situation as if we had merely classical igno- 
rance about the position of the particle. Needless to say, care 
is required in the interpretation of the result, and one should 
not be tempted by a realistic interpretation in terms of classi- 
cal alternatives. In turn, the total state of both the system and 
its environment is very different in structure and is typically a 
highly correlated and often, but not necessarily Q, entangled 
state. It is the hope that this paper can contribute to the debate 
on the dynamical appearance of classical properties in quan- 



tum theory. This debate is potentially becoming more timely 
than ever with the availability of novel experiments on deco- 
herence fl(ill . let it be with microwave cavities, ion traps, or 
nano-electromechanical systems. 

I would like to thank W.H. Zurek, H.-P. Breuer, M. Cramer, 
J.J. Halliwell and C. Henkel for very thoughtful and de- 
tailed comments on the manuscript and M.B. Plenio and 
P. Hanggi for discussions. This work was supported by 
the EU (QUPRODIS, QUIPROCONE), the A.v.-Humboldt- 
Foundation, the ESF, and the DFG. 



[1] W.H. Zurek, Rev. Mod. Phys. 75, 715 (2003); J.-P. Paz and 
W.H. Zurek, quant-ph/001001 1 D. Giulini, E. Joos, C. Kiefer, 
I.-O. Stamatescu, and H.D. Zeh, Decoherence and the Appear- 
ance of a Classical World in Quantum Theory (Springer, Hei- 
delberg, 1996); D. Guilini, C. Kiefer, and H.D. Zeh, Phys. Lett. 
A 199, 291 (1995); J.R. Anglin, J.P Paz, and W.H. Zurek, 
Phys. Rev. A 55, 4041 (1997); H.-P. Breuer and F. Petruccione, 
The Theory of Open Quantum Systems (Cambridge University 
Press, Cambridge, 2002). 

[2] A.O. Caldeira and A.J. Leggett, Physica A 121, 587 (1983). 

[3] N. Gisin and I.C. Percival, J. Phys. A 26, 2233 (1993); J.J. Hal- 
liwell and A. Zoupas, Phys. Rev. D 52, 7294 (1995). 

[4] D. Braun, F. Haake, and W.T. Strunz, Phys. Rev. Lett. 86, 2913 
(2001). 

[5] J. Eisert and M.B. Plenio, Phys. Rev. Lett. 89, 137902 (2002); 
H. McAneney, J. Lee, and M.S. Kim, Phys. Rev. A 68, 063814 
(2003). 

[6] W.H. Zurek, S. Habib, and J.P. Paz, Phys. Rev. Lett. 70, 1187 
(1993). 

[7] L. Diosi and C. Kiefer, J. Phys. A 35, 2675 (2002). 
[8] S. Bose, K. Jacobs, and P. L. Knight, Phys. Rev. A 59, 3204 
(1999); A.V. Rau, J.A. Dunningham, and K. Burnett, Science 
301, 1081 (2003). 
[9] F. Haake and R. Reibold, Phys. Rev. A 32, 2462 (1985); W.G. 
Unruh and W.H. Zurek, Phys. Rev. D 40, 1071 (1989); R. Kar- 
rlein and H. Grabert, Phys. Rev. E 55, 153 (1997). 

[10] S. Haroche, Physica Scripta T76, 159 (1998); Q.A. Turchette 
et al., Phys. Rev. A 62, 053807 (2000); DA. Kokorowski, A.D. 
Cronin, T.D. Roberts, and D.E. Pritchard, Phys. Rev. Lett. 86, 
2191 (2001); K. Hornberger, S. Uttenthaler, B. Brezger, L. 
Hackermiiller, M. Arndt, and A. Zeilinger, ibid. 90, 160401 
(2003); X.M.H. Huang, C.A. Zorman, M. Mehregany, and M.L. 
Roukes, Nature 421, 496 (2003). 

[11] B.L. Hu, J.P. Paz, and Y. Zhang, Phys. Rev. D 45, 2843 (1992). 

[12] V. Hakim and V. Ambegaokar, Phys. Rev. A 32, 423 (1995). 

[13] W. Weidlich and F. Haake, Z. Phys. 185, 30 (1965); P. Ullersma, 
Physica 32, 27 (1966); H. Grabert, U. Weiss, and P. Hanggi, 
Phys. Rev. Lett. 52, 2193 (1984); PS. Riseborough, P. Hanggi, 
and U. Weiss, Phys. Rev. A 31, 471 (1985); P. Talkner, Ann. 



Phys. 167, 390 (1986). 

[14] H. Grabert, P. Schramm, and G.L. Ingold, Phys. Rep. 168, 155 
(1988); H. Grabert, P. Schramm, and G.L. Ingold, Phys. Rev. 
Lett. 58, 1285 (1987). 

[15] Note that it is an interesting problem in its own right to investi- 
gate the influence of correlated initial conditions, such that time 
evolution of the state of the distinguished system no longer cor- 
responds to a completely positive map. This is relevant, e.g., 
when considering a charged quantum particle coupled to the ra- 
diation field, where the assumption of product initial conditions 
is typically not invariant under gauge transformations Hill . 

[16] H.-P. Breuer and F. Petruccione, Phys. Rev. A 63, 032102 
(2001). 

[17] GW. Ford and R.F. O'Connell, Phys. Rev. D 64,105020 (2001). 
[18] The convention for the Wigner function is 

with a being the symplectic matrix, such that £cr£' T is the stan- 
dard symplectic scalar product between . 

[19] For simplicity of the argument, the mass of the particle will set 
to unity and units will be employed where h = 1 and k = 1. 

[20] C. Sparber, J.A. Carrillo, J. Dolbeault, and PA. Markowich, On 
the long time behavior of the quantum Fokker-Planck equation, 
Monatshefte Mathematik 141, 237 (2004). 

[21] M. Hillery, R.F. O'Connell, M.O. Scully, and E.P Wigner, 
Phys. Rep. 106, 121 (1984). 

[22] This is meant in the sense that for any e > there exists a £ > 
such that |Boo -T| < e. 

[23] Needless to say, there is a certain degree of ambiguity in the 
choice of the pointer states here. The key argument is that the 
final state can after a finite time be exactly represented as a con- 
vex combination of these reasonable pointer states. 

[24] W. Schleich, Quantum Optics in Phase Space (Wiley- VCH, 
Berlin, 2001). 

[25] GM. D'Ariano and M.F. Sacchi, Nuovo Cim. B112, 881 
(1997). 



